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Abstract 

We study elastic proton-nucleus scattering at intermediate energies. The 
nucleon-nucleus optical potential is derived from the Bonn nucleon-nucleon 
potential and the Dirac-Brueckner approach for nuclear matter. Our calcula- 
tions, which do not contain any adjustable parameters, yield good agreement 
with experiment. 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



The microscopic description of the bulk properties of nuclear matter, finite nuclei and 
nuclear reactions in terms of a realistic nucleon-nucleon (NN) interaction continues to be an 
interesting topic in nuclear physics. There are two aspects to the problem. First, one needs a 
realistic NN potential. This strong interaction should, in principle, be derived from quantum 
chromodynamics (QCD). However, up till now, attempts to derive the NN interaction from 
QCD have led only to very crude agreement with experiment which is unsuitable for studies 
of nuclear many-body systems. Presently, the only quantitative models for the nuclear force 
are based on meson exchange. Well-known examples for such realistic potentials are the 
Bonn and the Paris fU potential; we will apply the former in this work. 

The second aspect of the problem is a suitable many-body approach that is able to deal 
with a realistic NN interaction which has strong short-range components. The Brueckner 
approach and the variational method have been developed for this purpose. However, when 
using two-body forces in these non-relativistic approaches, one is faced with a fundamental 
problem: even the bulk properties of nuclear matter cannot be reproduced correctly @,§]]. 

Encouraged by the success of the Dirac phenomenology for nucleon-nucleus scattering 
IIU and the Walecka model for dense nuclear matter [7||§, a relativistic extension of the 
Brueckner method has been initiated by Shakin and co-workers P,|TU| , frequently called the 
Dirac-Brueckner-Hartree-Fock (DBHF) approach. This approach has been further developed 
by Brockmann and Machleidt [|TTl|l2f and by ter Haar and Malfliet fl3| . Horowitz and Serot 
have discussed in detail the basic aspects involved in the derivation of the relativistic G- 
matrix |L4],[15|. The common feature of all DBHF results is that a repulsive relativistic 
many-body effect is obtained which is strongly density dependent such that the empirical 
nuclear matter properties can be explained, starting from a realistic NN interaction. 

There are basically two motivations for the extension of the DBHF approach to other 
observables of nuclear physics, e.g., nucleon-nucleus scattering. The first motivation is fun- 
damental: one wants to test if it is possible to describe nuclear matter, finite nuclei and 
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nuclear reactions in terms of the same basic NN interaction. Secondly, the relativistic de- 
scription of intermediate-energy nucleon-nucleus scattering is by itself an interesting topic 
of theoretical nuclear physics. Its most characteristic feature is the appearance of a strong 
attractive scalar and a strong repulsive vector potential, in terms of which the single-particle 
motion of the projectile nucleon in the mean field of the target nucleus can be described 
using the Dirac equation. In the Dirac phenomenology developed by Clark and collabo- 
rators H,|l^,|l7[] , these potentials are adjusted to the nucleon-nucleus scattering data, using 



an appropriate number of parameters. Besides the Dirac phenomenology, approaches have 
been developed in which nucleon-nucleus scattering is described more microscopically. We 



mention here the relativistic impulse approximation [p!q -pO| and calculations based on the 
Walecka model 

These approaches have some disadvantages. The Walecka model has no connection to 
the free-space NN interaction. In the relativistic impulse approximation, the medium effects 
are ingored. Thus it is tempting to develop a method that avoids both of these drawbacks. 
The DBHF calculation based on a realistic NN interactions offers this opportunity. There 



has been some work in this direction, e.g., in Refs. [0,0] the DBHF equations are solved 
directly in nuclear matter for a nucleon above the Fermi level. 

In this paper we present an alternative approach which combines the microscopic origin 
of the DBHF calculation with the simplicity of the Walecka model, such that an investigation 
of intermediate-energy nucleon-nucleus scattering is feasible. 

The direct extension of the DBHF calculation from nuclear matter to finite nuclei and 
nuclear reactions is not trivial. More practically, one can parameterizes the DBHF results 
for nuclear matter in terms of an effective Lagrangian which leads to the same predictions 
for nuclear matter as the original DBHF calculation. The effective Lagrangian, with its 
parameters derived from the realistic NN interaction, can then be used to investigate the 
properties of finite nuclei and nuclear reactions. Different schemes for this parameterization 
have been proposed We use in the present work a scheme similar to the one which 



was recently suggested by Brockmann and Toki ||28|| : The DBHF results for nuclear matter 



obtained from the Bonn potential are parameterized in terms of an effective a-u Lagrangian 
allowing for density-dependent meson-nucleon coupling constants. 

We discuss in Section II the details of our calculation. The results and a discussion are 
presented in Section III, Finally a short summary is given in Section IV. 



II. DETAILS OF THE CALCULATION 



The one-boson-exchange (OBE) Bonn potential p|,^9f used in this work is constructed 



in terms of the Thompson equation. The kernel of this integral equation is the sum 
of one-meson-exchange amplitudes of six nonstrange bosons with given mass and coupling. 
Pseudovector (derivative) coupling is used for pseudoscalar mesons (tt and rf) . A form factor 
of monopole type is applied to each meson-nucleon vertex which simulates the short-range 
dynamics of quark-gluon nature. See Ref. |2|j29|j for details about the OBE Bonn potential 
and its description of the two-nucleon system. 

As in conventional Brueckner theory, the basic quantity in the DBHF approach is a 
reaction matrix, G, which satisfies the in-medium Thompson equation. The nuclear matter 
properties are then obtained from this effective in-medium two-body interaction. Using the 
Bonn A potential ||, the DBHF calculation predicts that nuclear matter saturates at a 
density p=0.185 fm -3 and an energy per nucleon of -15.6 MeV, which is in good agreement 
with empirical information. More results and discussion concerning the properties of nuclear 
matter as obtained in the DBHF approach with the Bonn potential have been presented in 
Refs. |2|,|l2|,[n|] . In this work, we use the Bonn A potential 



As proposed by Brockmann and Toki [28], the DBHF results for nuclear matter can be 



parameterized by an effective Lagrangian, in analogy to the o-uj model of Walecka: 
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L = -m- g a (p)<p a - gM^M + -(d^ a ) 

- ~(<^C - + (i) 
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where ip is the nucleoli field, while CT and 0^ are effective sigma and omega fields, respec- 
tively. The masses of the sigma and omega mesons, m a and m w , resp., are kept fixed at their 
values in free-space scattering (m a = 550 MeV, m w = 782.6 MeV). However, the coupling 
constants of these effective mesons, g a and g^, depend on the baryon density p. They are 
determined from the DBHF results for nuclear matter. 

The (complex) nucleon self-energy produced by the effective meson exchanges in nuclear 
matter can, in general, be written as 



where £5, So and denote the scalar component, the time-like part of the vector compo- 
nent and the space-like part of the vector component of the nucleon self-energy, respectively. 

In our notation, k^ stands for all four components of the momentum four-vector, and 
k = |k|. 

Treating the effective coupling constants locally as numbers and calculating in the rel- 
ativistic Hartree-Fock approximation (see Fig. la and lb for the corresponding Feynman 
diagrams), we obtain the explicit expressions for the real part of the nucleon self-energy 



Vs(k, p) = -i^fPs + dqq^j^tfMe^k, q) - 4^(p)0^, q)} (3) 

Vb(M = (iM)2 p+z J_ f F d qq {gl{p)Q a {k, q ) + 2gl(p)Q w (k, q )} (4) 
m u lo7T K Jo 



£(Ag = S s (^) + 7°So(M + 7 ■ kE v (^) =V + iW 



(2) 



mum 



V v {k,p) 



8n 2 k 2 Jo 




where 




(27r)»£lj(p) 
r k F d 3 q 




and 



Qi(k,q) = ln\ 



(k + qf + m 2 - (q - k ) 2 
(k - q) 2 + m 2 - (g - k ) 2 
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*n \ k 2 + q 2 + mf - (k - q ) 2 „ n , . 
$(fc,g) = — Qi{k,q)-1, i = a,u 

The effective mass m*, the kinetic momentum k* and the single particle energy ko are given 
by 

m l(p) =m + V s (k,p) 
k* = k(l + V v (k)) 
k = El + V = (k* 2 + mf ) 1/2 + V 

and similarly for q. 

The first term of V$ and Vq is the momentum-independent Hartree contribution (see Fig. 
la), while the other terms are the Fock contributions (see Fig. lb) which depend on k. A 
very small space-like component of the vector potential is due entirely to Fock diagrams. 

In order to determine the effective coupling constants from these expressions, we drop the 
space-like component of the vector potential (Eq. (5)) and calculate Vs (Eq. (3)) and Vq (Eq. 
(4)) for k = kp, which are then identified with the scalar and vector potentials obtained in 



the DBHF calculation (see Table VII of Ref. |12[)- This is a reasonable assumption since the 



space-like component of the vector potential is very small and the potentials are only very 



weakly momentum dependent ||33|| . Note that these approximations are used only for the 
determination of the effective coupling constants; in the calculation of the optical potential, 
the space-like part of the vector potential (Eq. (5)) is taken into account (see Eq. (10) 
below). The momentum k in Eqs. (3)- (5) is related to the incident energy of the projectile 
(see Eq. (11) below) through ko = E. Both effective coupling constants drop with increasing 
density 

The lowest-order contribution to the imaginary part of the nucleon self-energy is obtained 
from the fourth-order (counting the meson-nucleon vertices) Feynman diagrams which are 
characterized by two-particle-one-hole (2plh) intermediate states (see Fig. lc). The nucleon 
lines in these Feynman diagrams are described by dressed nucleon propagators, which cor- 
responds to performing the calculation on the Hartree-Fock ground state. The imaginary 
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part of nucleoli self-energy due to the fourth-order Feynman diagrams can be expressed as 



W s(k,p) = — 7pr / dq dqq[g 2 a {p)m*Al{q)Imir a {q) 
6n z k Jo J\k-p\ 



2 

+g 2 Mm*Al(q) gflX lmir^(q) + 2 9er (p)gM(E^ - \q )%K{q)K{q)l^lM\ (6) 



I f E' f k+p 

^o(M = 77-27;/ dqj dqq\g 2 MAl(q){EZ-q )Tmir„(q) 

8-K z K JO J\k-v\ 



8ir 2 k Jo J\k-p\ 

2 

-,2 J 



-9 2 M^l(Q)m - qo)9,xlmn^(q) + (2E* k - qo)^!™™ (q)} 



-2g a (p)g^p)m*A a (q)AUq)Im7rL(q)} (?) 

1 pE 1 rk+p I 

W v(k,p) = -tt^u / dq dqq[gl(p)Al(q)(k 2 + -q 2 - E^q ))lm7T a (q) 

5H*k A Jo J\k-p\ 2 p 

+gUp)^l(q)i-(k 2 + \i - ^o)^Almvr^(g) + ^(2£^ - q 2 ^ - E* h )lnvK™ (q)} 

-g a (p)g^ p ) m *^(2k*q - fi)A o (q)A u {q)lm<n au {q)] (8) 
where q is the magnitude of the three-vector q, while 

m * = m l(p) 
E' = E*-(k 2 F + m* 2 ) 1/2 
p ={{Et-q,) 2 -m* 2 f/ 2 

Ql = %(f 
A t (q) = -=— — j, % = a, uj 

% - m 

7r(q) are meson polarization insertions; their imaginary parts are given below 

Imvr CT (g) = - 9 -M(m* 2 - l^)(ET F - E*) 
irq 4 ^ 

Imvr° (g) = - 9 ^(\{Ef - E* 3 ) + \{Ef - E* 2 )q + ]q 2 (E F - E*)) 
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WL(g) = - 9a{p } 9M m*{{Ef - £* 2 ) + (E F - E*)q ) 



where 



E* = max{E r , E* F — q , m*} 
E r = ~\qo + \q{l ~ ^/qlf/\ E F = (k 2 F + m* 2 ) 1 ^ 
The derivation of the nucleon self-energy (Eqs. (3)-(8)) from the Walecka model has 



been discussed in detail in Refs. |2lj22 l. For the effective Lagrangian (Eq. (1)) used in 
the present work, the expressions for the nucleon self-energy are the same, only with the 
coupling constants of the sigma- and omega-exchanges replaced by those determined from 
the DBHF approach which are now density dependent. 

The same fourth-order diagrams, which yield the imaginary part also give rise to a 
contribution to the real part of the self-energy, which should be added to the Hartree- 
Fock term. Note, however, that the Hartree-Fock contribution has been derived from a 
G-matrix description of nuclear matter and therefore already contains these terms in a 
certain approximation. In order to avoid any double counting no contributions of the fourth 
order diagram to the real part have been considered. 

The nucleon self-energy in a finite nucleus is obtained by means of the local density ap- 
proximation, in which the spatial dependence of the microscopic optical potential is directly 
related to the density of the target nucleus. For a self-consistent description of nucleon- 
nucleus scattering (by self-consistency we mean that there are no free parameters after the 
NN potential has been selected), the target density should also be determined from the ef- 
fective Lagrangian, Eq. (1). In this paper, we use the density determined in the relativistic 
density-dependent Hartree-Fock (RDHF) calculation with the effective Lagrangian of Eq. (1) 



p5| . For comparison, we also use the density determined in a relativistic density-dependent 



Hartree (RDH) calculation |35[, as well as the density determined in the relativistic Hartree 



(RH) calculation with the Walecka model |36fl . As shown in Ref. |[35|| , the RDH results 
for the binding energy and root-mean-square radius of 40 Ca are in good agreement with 
experiment, while the RDHF results slightly underestimates them. 
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The Dirac equation for the single-particle motion of the projectile nucleon in the mean 
field of the target nucleus can be written as 



[a-k + p{m + U s ) +U V + V C 



c 



with 



Uc 



U 



S + 



V 



(9) 



(10) 



where E is the energy of the projectile in the center-of-mass (cm.) system of projectile and 
nucleus, which is related to the incident energy 7] a & through 

m 2 + m-rfm + TJ C 



E = E. 



-lab 



p,c.m. 



(11) 



((m + m T y + 2m T T lab y/ 2 
with m and the mass of the projectile and target, respectively. Vc is the Coulomb field 
which we treat in the same way as in Refs. [2(J,36 . 

In order to calculate the experimental observables, the Dirac equation is usually converted 
into a Schrodinger-equivalent equation. The Dirac equation for the four-component spinor if) 
is equivalent to two coupled equations for the large (upper) and small (lower) two-component 
spinors. One can eliminate the small component of the Dirac spinor in a standard way 
and obtain a Schrodinger-equivalent equation for the large component of the Dirac spinor 



2E 



+ U eff {r) + V c {r) + U s . .{r)a-1 



E 2 



m 



2E 



(12) 



where C/ e // and U s , , are the central and the spin-orbit part of the Schrodinger-equivalent 
potential, which is known as the nucleon-nucleus optical potential in non-relativistic ap- 
proaches. The explicit expressions for U e ff and U s _ _ are 



U eff = U V + ^[U s (2m + U s ) - (U v + V c f + U Darwm ] 

1 dD(r) 



2E l 

U s .o. 



2ErD(r) dr 



(13) 
(14) 



with 



U 



Darwin 



1 dD(r) 
D(r) dr 



1 2 



dD(r) 



d 2 D(r) 



rD(r) dr 



2D(r) dr 2 



D(r) = m + E + Us-U v -V c 



c 



III. RESULTS AND DISCUSSION 



We show in Fig. 2 the density of 40 Ca determined in a RDHF (solid curve) and RDH 
(dashed curve) calculation using the effective Lagrangian of Eq. (1); they are compared 
with the density (dotted curve) as obtained from a relativistic Hartree (RH) calculation 
applying the Walecka model (QHD-I) |ftfl|| . Note that the RDHF calculation slightly under- 
estimates the root-mean-square radius of 40 Ca PJ. 

We show in Fig. 3 the scalar potential Us and vector potential Uy for p+ Ca scattering 
obtained in the present calculation. The RDHF density is used in the local density approxi- 
mation. The real parts of these potentials are only weakly energy dependent; in the interior 
of the nuclei, they are about —420 and +320 MeV, respectively. The overwhelming part of 
these real potentials comes from the energy independent Hartree contribution; the energy 
dependent Fock contribution plays only a minor role. On the other hand the imaginary 
part of the scalar and vector potential depends of course very strongly on the energy; both 
increase in magnitudes with increasing energy. Note that the radial shapes of the real Dirac 
potentials follow essentially that of the target density (cf. solid curve in Fig. 2). 

With the nucleon self-energy determined from the effective Lagrangian which is related 
to the realistic NN interaction used in the DBHF calculation, we can obtain the Schrodinger- 
equivalent potential through Eqs. (13) and (14) which is then used in Eq. (12) to calculate 
the experimental observables f36|. The Schrodinger-equivalent potentials, corresponding to 
the Dirac potentials in Fig. 3, are shown in Fig. 4. As can be seen from the figure, the real 
part of the central potential is strongly energy dependent. At low incident energy (Ti a f,=150) 
the real part of the central potential is attractive throughout the whole radial space and has 
a longer range. With the increase of the incident energy, the real part of the central potential 
becomes less attractive and is repulsive in the interior of the nucleus when T; a 6=300 MeV. 
A typical pocket of attraction is observed at the nuclear surface at this energy. At even 
higher energy (7] a &=450 MeV), the real central potential is repulsive in the whole radial 
space. The imaginary part of the central potential also shows a strong energy dependence. 
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The imaginary central potential is negative at all energies; its magnitude increases while its 
range decreases with increasing energy. 

The important feature of the relativistic approach is that the spin-orbit potential arises 
naturally from the coherent sum of the contributions from the scalar and vector potential. 
While the real part of the spin-orbit potential is attractive, its imaginary part is mostly 
positive. The real part of the spin-orbit potential decreases (in magnitude) with increasing 
incident energy, whereas its imaginary part increases with the energy. The fluctuations of 
the spin-orbit potential in the interior of the nucleus are due to fluctuations in the target 
density (cf. solid curve in Fig. 2). 

The results for differential cross sections and analysing powers in elastic p+ 40 Ca scat- 
tering at T/ a f,=300 and 400 MeV are shown in Fig. 5 and Fig. 6, respectively. In these 
figures, the solid and dashed curves are the results of our calculations with RDHF and RDH 
densities, respectively, while for the dotted curves the RH density from the Walecka model 
(QHD-I) |36[] is used. The experimental data are from Ref. ||37|| . 

For full self-consistency of our calculations, we must use the RDHF density since it 
is based upon the same effective coupling constants as our optical potential. With the 
RDHF density, our results (solid curves) are in reasonable agreement with experiment. The 
agreement is particularly good at small angles (8 c , m , < 15°). At large angles, we overestimate 
the experimental differential cross section. The oscillations in the differential cross section 
and analysing power seem to have a larger period than the experimental data show (note that 
there is a correspondence between the minima in the differential cross section and the dips 
in the analysing power). Since our calculations are parameter free, the quality of agreement 
with the experimental data is remarkable. The remaining discrepancies are probably mainly 
due to the fact that the target density as obtained in the RDHF calculation underestimates 
the root-mean-square radius J£J (see discussion below). 

When using the RDH density, which is in better agreement with the experimental root- 
mean-square radius our results (dashed curves in Figs. 5 and 6) improve; and there is 



further improvement when the RH density (dotted curves) is used. Thus, the accurate re- 
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production of the nucleon density of finite nuclei is probably the most important outstanding 
problem for the self-consistent description of nuclear structure and nuclear reactions. 

IV. SUMMARY 

In this paper, we have continued our efforts of obtaining a parameter-free and selfcon- 
sistent description of nuclear matter, finite nuclei, and nuclear reactions in terms of one 
realistic NN interaction. Here, we have studied p+ 40 Ca scattering in Dirac dynamics. In 
order to carry out a systematic study, we parameterized the DBHF results for nuclear mat- 
ter in terms of a simple effective Lagrangian consisting of an attractive sigma-exchange and 
a repulsive omega-exchange. The coupling constants of these effective mesons are density 
dependent and are determined from the DBHF results for nuclear matter. 

The nucleon self-energy in nuclear matter is then calculated based on the effective La- 
grangian up to the fourth order Feynman diagrams. The nucleon-nucleus optical potential in 
finite nuclei is related to the self-energy in nuclear matter by means of the local density ap- 
proximation, where the target density is determined by the same effective Lagrangian as the 
nucleon self-energy. Furthermore, the Dirac equation describing the single-particle motion 
of the projectile in the mean field of the target is converted to the Schrodinger equivalent 
equation which is solved to obtain the observables. 

With the optical potential thus determined, we have calculated observables of p+ 40 Ca 
scattering at intermediate energies. The predictions by our model for the elastic differential 
cross section and the analyzing power were compared with the experimental data. The 
agreement between our parameter-free calculations and experiment is remarkable. At large 
scattering angles, our results for the differential cross section overestimate the experimental 
data. This can be attributed, in part, to small deficiencies in the densities as obtained in 
the RDHF calculation. 
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Figure Captions 



Fig. 1: Feynman diagrams for the calculation of the nucleon self-energy in nuclear 
matter, (a) Hartree contribution , (b) Fock contribution , (c) fourth-order contribution. 

Fig. 2: Nucleon density of 40 Ca as obtained in RDHF (solid curve) and RDH (dashed 
curve) calculations with the effective Lagrangian, Eq. (1) P5|; the dotted curve represents 
a RH calculation with the Walecka (QHD-I) model 

Fig. 3: Scalar and vector potentials in p+ 40 Ca scattering at 150 (solid curves), 300 
(dashed curves) and 450 MeV (dotted curves). The RDHF density is used in local density 
approximation. 

Fig. 4: The Schrodinger-equivalent potential in p+ 40 Ca scattering at 150 (solid curves), 
300 (dashed curves) and 450 MeV (dotted curves). The RDHF density is used in local 
density approximation. 

Fig. 5: Elastic differential cross section and analyzing power in p+ 40 Ca scattering at 
21a6=300 MeV. The solid and dashed curves are the predictions by the present calculations 
with RDHF (solid curve in Fig. 2) and RDH (dashed curve in Fig. 2) densities, respectively, 
while the dotted curves are the results using the RH density from the Walecka (QHD-I) 



model (dotted curve in Fig. 2). The experimental data are from Ref. |37 
Fig. 6: Same as Fig. 3, but at T lab =AQQ MeV. 
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